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A microscopic theory of the phonon-induced resistance oscillations in weak perpendicular magnetic 
fields is presented. The calculations are based on the consideration of interaction of two-dimensional 
electrons with three-dimensional (bulk) acoustic phonons and take into account anisotropy of the 
phonon spectrum in cubic crystals. The magnetoresistance is calculated for [001]-grown GaAs 
quantum wells. The results are in agreement with available experimental data. Apart from the 
numerical results, analytical expressions for the oscillating part of magnetoresistance are obtained. 
These expressions are valid in the region of high-order magnetophonon resonances and describe 
the oscillating magnetoresistance determined by several groups of phonons polarized along certain 
high-symmetry directions. 
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I. INTRODUCTION 

In recent years, experimental studies of transport prop- 
erties of high-mobility two-dimensional (2D) electron gas 
in weak perpendicular magnetic fields have uncovered a 
variety of remarkable quantum phenomena caused by 
transitions of electrons between different Landau lev- 
els. Such transitions can lead to oscillations of dissipa- 
tive resistance as a function of the magnetic field. For 
example, under steady-state microwave illumination of 
the 2D gas the resistance oscillates with a period de- 
termined by the ratio of the radiation frequency to the 
cyclotron frequency lOc- This phenomenon is known as 
the microwave-induced resistance oscillations (MIRO).^ 
For a sufficiently high radiation power the MIRO min- 
ima evolve into the intervals of magnetic field where the 
dissipative resistance vanishes. ^''^ Next, it was found that 
an increase in the electric current passing though the 2D 
layer substantially reduces the resistance^'^ and leads to 
oscillations of the resistance as a function of either the 
magnetic field or the current. ^'^'^ Such oscillations are 
controlled by the ratio of a characteristic energy, which 
is defined as the drop of the Hall electric field across 
the classical cyclotron diameter, to the cyclotron energy. 
This phenomenon has been called the Hall field-induced 
resistance oscillations (HIRO). It is also described in 
terms of Zener tunneling between Landau levels. Finally, 
in the systems with two occupied 2D subbands the resis- 
tivity oscillates as a function of the ratio of subband sep- 
aration energy to the cyclotron energy owing to coupling 
of the subbands via intersubband scattering. This phe- 
nomenon, called the magneto-intersubband oscillations 
(MISO), is important in the systems with small subband 
separation, such as the double quantum wells. All these 
kinds of oscillations are insensitive to positions of the 
Landau levels with respect to the Fermi energy. There- 
fore, unlike the Shubnikov-de Haas oscillations, they are 
not damped exponentially with increasing temperature. 

The oscillatory phenomena described above require the 
presence of microwave illumination, strong Hall field, or 



intersubband coupling to enable resonant transitions be- 
tween different Landau levels via elastic scattering of elec- 
trons by impurities or other static imperfections of the 
system. The elastic scattering provides the main con- 
tribution to resistivity at low temperatures. In high- 
mobility samples, however, inelastic scattering of elec- 
trons by acoustic phonons also contributes to trans- 
port, and can exceed the impurity-scattering contri- 
bution at the temperatures of several K. Therefore, 
acoustic-phonon-assisted transitions of electrons between 
Landau levels become important. Interestingly enough, 
these transitions in 2D case possess a resonant property 
and, for this reason, also lead to oscillations of the resis- 
tance. Oscillations of this origin have been discovered 
recently. An analysis of experimental data has sug- 
gested that the resistivity acquires an oscillating contri- 
bution whose periodicity is determined by the resonant 
condition 

2pps — nhuJc, (1) 

where pp is the Fermi momentum of electrons, s is the 
constant of the order of sound velocity, and n is an 
integer. This phenomenon has been called the phonon- 
induced resistance oscillations (PIRO). Similar oscilla- 
tions have been observed later in the phonon drag ther- 
mal power measurements.^^ 

Since the quantity in the left-hand side of Eq. (1) is 
a characteristic phonon energy, the PIRO can be viewed 
as magnetophonon oscillations. For the case of electron 
scattering by optical phonons, the magnetophonon os- 
cillations are known for a long time^^ and have been 
observed both in bulk and 2D-'^'^ systems. The special 
feature of PIRO, in contrast to the oscillations owing 
to scattering by optical phonons, is the dependence of 
magnetophonon resonance on the electron momentum. 
By emitting or absorbing acoustic phonons, the elec- 
trons jump between the Landau levels, and the max- 
imum probability of such transitions is realized under 
electron backscattering condition, when the phonon wave 
number q reaches its maximum value 2pp/h in the 2D 



2 



plane. Initially, -'^'^ this resonant property of the electron- 
phonon scattering was explained by involving the inter- 
face phonon model. Later, by analyzing phonon-assisted 
transitions between Landau levels, it was found that 
the bulk phonon model gives similar results. It is worth 
noting that the maximum probability of the backscatter- 
ing processes {q = 2pp/h) is not related to the pres- 
ence of the magnetic field. It reflects a fundamental 
property of the kinematics of two-dimensional electron 
scattering by thre e-dimen sional acoustic modes with fre- 
quency Wq = s^J q]_ + q^, where q±_ and qz are the in- 
plane and out-of-plane components of the phonon wave 
vector. It can be shown that the scattering probability, 
as a function of the phonon frequency, has a logarith- 
mic singularity at the point of transition from the region 
Wq < 2pps/h, when pure 2D scattering [q^ = 0) is possi- 
ble, to the region Wq > 2pFs/h, when emission of three- 
dimensional phonons (with finite Qz) becomes necessary. 

The systematic studies of the PIRO are now un- 
der way.^^"-^'^ However, current understanding of this 
phenomenon is far from complete. The characteristic 
phonon velocity s entering Eq. (1) varies in different 
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in the range from 2.9 to 5.9 km/s. 



and the origin of such variations is not clear. The prob- 
lem of the phase of the oscillations has not been dis- 
cussed. The dependence of the amplitude of the oscil- 
lations on temperature and magnetic field has not been 
investigated in detail, though general features of the ex- 
perimental dependence have been successfully related"'^^ 
to the behavior of the density of electron states in mag- 
netic field. The main difficulty for interpretation of the 
experimental magnetoresistance is the complicated na- 
ture of acoustic-phonon modes. Even in the case when 
the influence of interfaces on the phonon spectrum is ne- 
glected (bulk phonon approximation), one should take 
into account three different phonon branches character- 
ized by anisotropic (i.e., dependent on the direction of 
phonon momentum) velocities. 

In this paper, the theory of PIRO is developed within 
the bulk phonon approximation. Both deformation- 
potential and piezoelectric mechanisms of electron- 
phonon interaction are included into consideration. The 
anisotropy of phonon spectrum is taken into account. 
The calculations are carried out for GaAs quantum 
well layers grown in the [001] crystallographic direction. 
Apart from the numerical results, approximate analytical 
expressions for the oscillating part of the phonon-induced 
resistivity are presented. The results are compared with 
available experimental data for two kinds of GaAs lay- 
ers with distinctly different densities and mobilities of 2D 
electrons. 

The paper is organized as follows. Section II contains 
general description of the resistivity of 2D electrons in- 
teracting with acoustic phonons in a weak magnetic field. 
Section III includes both analytical and numerical results 
for the resistivity, their comparison with experiment, and 
discussion. The summary and conclusions are given in 
the last section. 



II. GENERAL FORMALISM 

Experimental investigations of phonon-induced resis- 
tance oscillations are carried out at low temperatutes 
T satisfying the condition of degenerate electron gas, 
T <^ ep, where ep is the Fermi energy. Magnetic fields 
are weak enough to have many Landau levels populated, 
hujc <^ ep, so the electron transport possesses some qua- 
siclassical features (for example, the electron momentum 
is still a good quantum number). In the same time, 
the high quality of the samples ensures that the Lan- 
dau quantization is not suppressed by scattering at such 
low magnetic fields. This means that the ratio of the cy- 
clotron frequency ujc to inverse quantum lifetime 1 /r^ is 
not too small. Next, since in the high-mobility samples 
the ratio of transport time to quantum lifetime is large, 
the oscillations are observed in the regime of classically 
strong magnetic fields, when the cyclotron frequency is 
much larger than the inverse transport time l/rtr- The 
in-plane transport in this case can be viewed as hopping 
between cyclotron orbit centers, and the dissipative con- 
ductivity ad is represented in the form 

'^ci = E/(^^)[i - fMm'iXs - x,,fi2, (2) 
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where (5 is the index of a quantum state of electron, 
is the coordinate of the cyclotron orbit center in the di- 
rection of motion, /(e) is the equilibrium Fermi-Dirac 
distribution function, I? is the normalization area, and 
V is the total scattering rate, which is represented as 
a sum of impurity and phonon scattering contributions: 
= ^^55? + vf^i ■ Assuming interaction of electrons with 
equilibrium phonons, the latter contribution is written as 
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x^(e5 - £5/ - hjjJx^ + -/VAq5(£5 - £5' + ^\^\ , (3) 

where A, q = {qx,(iy,Qz), and cjxq are the phonon mode 
index, wave vector, and frequency, respectively. Next, 
-^Aq = [exp(fiwAq/T) — is the Planck distribu- 

tion function. The interaction is described by the func- 
tion Cxq, which includes both deformation-potential and 
piezoelectric coupling of electrons to crystal vibrations. 
Considering only acoustic modes in cubic crystals, the 
general form of this function is written as follows: 
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where V is the deformation potential constant, hu is 
the piezoelectric coupling constant, and p is the material 

density. The sums are taken over Cartesian coordinate 
indices, CAqi are the components of the unit vector of 
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the mode polarization, and the coefhcient Kijk is equal 
to unity if all the indices i,j,k are different and equal 
to zero otherwise. The polarization vectors and the cor- 
responding phonon mode frequencies are found from the 
eigenstate problem 

Yl [^ii ~ ^aP^"^] ^Aqj- = 0. (5) 

i 

The dynamical matrix for cubic crystals, 

+ (ci2 + C44)qtqj{l - (6) 

is written in the elastic approximation and expressed 
through three elastic constants for which the conven- 
tional notations are used. Equation (5) has three so- 
lutions describing a high-energy mode (LA), which be- 
comes purely longitudinal for high-symmetry directions 
([100], [110], and [111]), and two low-energy modes (TA) 
which become purely transverse for these directions. 
Very often, theoretical studies of transport properties 
of electrons in solids with cubic lattice symmetry are 
based on the continuum approximation for elastic vibra- 
tions, when C44 — (cii — ci2)/2 and the phonon spec- 
trum is given by three isotropic branches including one 
longitudinal mode and two degenerate transverse modes. 
Whereas this approach gives a reasonably good descrip- 
tion of phonon-limited mobilities and applies for calcu- 
lation of many kinetic coefficients, is not sufficient for 
description of phonon-induced resistance oscillations in 
magnetic field, because the anisotropy of the acoustic 
phonon branches is not weak in most semiconductors and 
becomes essential in evaluation of resonant scattering be- 
tween Landau levels. 

The quantum states 5 in Eq. (2) can be treated as 
exact eigenstates of electrons interacting with randomly 
distributed impurities in the presence of a magnetic field. 
It is convenient to rewrite Eq. (2) through the Green's 
functions of electrons. Using the basis of Landau eigen- 
states and specifying the growth axis of the quantum- well 
layer as z (i.e., the [001] crystallographic direction), one 
obtains the phonon-induced contribution to the conduc- 
tivity: 

X r|^E^Aq(A^Aq + l)E*-'(9i«l/2) 

X J defis)[l - f{e - hiJx^)]A,{n')A,_n^^^{n), (7) 

where q± and are the absolute value and polar angle of 

the in-plane component of phonon wave vector, n arc the 
Landau level numbers, and Ib is the magnetic length. 
The function ^nn'{x) = (n!/n'!)a;"'-"e-^[L;^'-"(a;)]2, 
where L'^{x) are the Laguerre polynomials, describes 
scattering in the magnetic field. The squared matrix 



element of a plane-wave factor, I (Ole^'-^IO) p, is deter- 
mined by the confinement potential which defines the 
ground state of 2D electrons, |0). Applying the model of 
a deep rectangular quantum well of width d^, one can 
rewrite this squared matrix element as I{qzdyj/2), where 
I{x) = (sin.x/x)^[l — (x/tt)^]"^. Finally, the spectral 
function As{n), which is equal to the imaginary part of 
the single-electron (advanced) Green's function divided 
by TT, characterizes disorder-induced broadening of elec- 
tron states. In the Born approximation, this broadening 
is conveniently described in terms of the quantum life- 
time of electron, Tq. In the absence of impurities (the 
coUisionless limit) (n) is reduced to the delta- function 
S{e — En), where e„ = hcodn + 1/2) is the Landau quan- 
tization energy. 

The dissipative resistivity measured in experiments us- 
ing Hall bars is related to the conductivity ad according 
to pd — ad/'^H^ where an = e^Us/moJc is the classi- 
cal Hall conductivity and Ug is the sheet carrier density. 
The resistivity is represented as a sum of contributions 
from electron-impurity and electron-phonon scattering: 

muim mvph 

Pd — Pirn + Pph, Pirn — — , Pph — — 5 \0) 

where Uim and Vph are the partial scattering rates which 
depend on the magnetic field. Both pph and Vph are 
straightforwardly obtained from of Eq. (7). If 
the magnetic field is zero, and are reduced to 
the inverse transport times due to electron-impurity and 
electron-phonon interactions, respectively, and Eq. (8) 
describes the classical (Drude) resistivity. 

Let us consider the case of weak magnetic fields, when 
many Landau levels are occupied. Then one can use 
smallness of the cyclotron energy fvjjc, phonon energy 
hiUXq, disorder- broadening energy h/Tq, and temperature 
T with respect to the Fermi energy for evaluation of the 
expression (7). Owing to the sharp form of the spec- 
tral functions, the contribution into the n-sums in Eq. 
(7) comes from a limited number of Landau levels in the 
vicinity of the Fcarni energy, and these simis can be ap- 
proximately calculated with the following result: 

2m /■^'^ dip /■^'^ d9 , 

X J def{s)[l - f{e - hioxq)]D{e)D{e - hvx^), (9) 

where D(s) is the dimensionless (i.e., normalized to its 
zero-magnetic-field value) density of electron states. The 
integration over the in-plane phonon wave number q± is 
replaced in Eq. (9) with the integration over the scat- 
tering angle according to q± = 2{pp/h)sm{0/2). This 
emphasizes that the phonon-assisted transport at large 
filling factors is described within the picture of quasielas- 
tic scattering of 2D electrons in the vicinity of the Fermi 
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surface. Under the same assumptions, the contribution 
to the resistivity owing to electron-impurity scattering is 
given by 



(10) 



where r/™ is the transport time for electron-impurity 
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scattering. It is defined according to l/r^, 

(m/27r?l^) dd w{2pF sin |)(1 — cos6), where w is the 
Fourier transform of the random impurity-potential cor- 
relator. The self-consistent Born approximation leads to 
the density of states in the form of an expansion in os- 
cillation harmonics weighted with powers of the Dingle 
factor d = exp{—TT /uJcTq):^^'^^ 
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flfc = (-1)''^ ^ exp(-7rA;/wcTg)ife_i(27rfc/wcTg). 



If the magnetic fields is weak, the Dingle factor is small. 
This corresponds to the case of overlapping Landau lev- 
els, when the density of states is given by a single- 
harmonic expression, D{e) = 1 — 2dcos{2TTe/fiu>c)- In 
strong enough magnetic fields the electron system is in 
the regime of separated Landau levels, when D{e) is rep- 
resented by a periodic sequence of semi-elliptic peaks cen- 
tered at the Landau quantization energies e„ . 

The expressions (8)-(ll), together with Eqs. (4)-(6) 
defining the phonon spectrum LUxq and the function C'xq, 
give a complete description of the phonon-assisted 2D 
magnetotransport at large filling factors. The next step 
in evaluation of the resistivity can be done by calculat- 
ing the integral over the energy e. In this procedure, 
it is reasonable to omit the terms containing oscillating 
functions of energy under the integral. Such terms arc 
responsible for the Shubnikov-de Haas oscillations, and 
they are exponentially suppressed with increasing tem- 
perature. Thus, the approximation used below corre- 
sponds to the condition 
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which means that the temperature should not be too low. 
As a result, 
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where F{x) = [a;/sinh(a;)]^ and Uz = hqz/2pF- The 
quantity 



hPojxq 



(14) 



is an anisotropic scattering rate introduced for con- 
venience purpose. The corresponding result for the 
impurity-assisted magnetotransport has a simpler form, 
^im = [1 + 2X]^i a^]/r(™, which describes a positive 
magnetoresistance oving to Landau quantization.^^ The 
Fermi momentum entering Eqs. (13) and (14) is related 
to electron density according to pp = h^/2Tml. 



III. RESULTS AND DISCUSSION 



The expression (13) is the central part of this paper. 
The first term of this expression, t'p^^ is responsible for 
phonon-induced contribution to background resistivity, 
while the second one corresponds to oscillations of the 
resistivity as a function of the magnetic field. These os- 
cillations (PIRO) are described in Eq. (13) through a 
sum of oscillatory harmonics of the scattering rate, i^^^^ . 
Because of the complicated form of phonon spectrum, the 
integrals in Eq. (13) cannot be calculated analytically in 
the general case. However, an analytical consideration is 
possible for the oscillating part of Vph, since this part con- 
tains rapidly varying functions of w^q under the integral. 
To carry out the approximate integration, it is conve- 
nient to represent expressions for oscillatory harmonics 
by using the spherical coordinate system according to 
sin(0/2) = usinx, Uz = ucosx, where u = hq/2pp is 
the absolute value of phonon wave vector in the dimen- 
sionless form and x is the inclination angle. In these 
variables, the phonon spectrum is conveniently written 
as uj\q = s'ip^q, where is the anisotropic sound ve- 
locity for the mode A. The angles ip and x determine the 
direction of the phonon wave vector. As a result. 
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where v\q given by Eq. (14) is written as the expression 
in the square brackets; the deformation-potential (DA) 
and piezoelectric (PA) contributions are separated ac- 
cording to Eq. (4). The main contribution to the inte- 
gral over u comes from the region u ~ (sinx)"'^, and the 
convergence of this integral is determined by the rapidly 
oscillating cosine function. The other w-dependent func- 
tions, / and _F, are slowly varying on the scale of conver- 
gence because oi ppdw/h ~ 1 and because of condition 
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(12), respectively. After integrating over u, one obtains 
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The remaining double integral over angular variables can 
be calculated using the method of fastest descent. In- 
deed, under the integral we have a rapidly oscillating 
function (the argument of the cosine in Eq. (16) is as- 
sumed to be large), and the main contribution comes 

from certain regions of ip and x where = ^vx I sin X 
varies most slowly. Naturally, these regions are in the 



close vicinity of the extrema points of the function 

The anisotropic sound velocity itself has a number 
of local maxima, minima, and saddle points, whose po- 
sitions coincide with some high-symmetry directions in 
the reciprocal (momentum) space. The number of ex- 
trema for >^\px is smaller because of the factor 1/ sinx. 
For the mode with the highest velocity (this mode be- 
comes purely longitudinal for high-symmetry directions) 
the anisotropy is not strong, and all the extrema are only 
at X = 7'"/2 (zero Qz). The slow modes (which become 
purely transverse for high-symmetry directions) also have 
extrema at x = 7r/2. Since the anisotropy of these modes 
is stronger, they also may have additional groups of ex- 
trema at X = 7i'/2 ± ^, where ^ < 7r/4. Nevertheless, the 
main contribution of slow (low-energy) modes to oscillat- 
ing magnetoresistancc is associated with the extremum 
X = 7r/2 for a single mode whose velocity at x = ''■/2 is 
independent of the polar angle ip. This mode at x = ""72 
is purely transverse and polarized perpendicular to the 
quantum well plane (direction [001]). 

In summary, the analysis shows that the modes ef- 
fectively contributing to the oscillating part of resistiv- 
ity are: i) the transverse mode polarized along [001], 
whose velocity Sto = ^Jc^aJp is independent of ip 
and corresponds to a maximum of <I>^^' as a function 
of x; ii) the longitudinal mode with velocity slo = 
ycujp polarized along [100] or along equivalent direc- 
tions ([010], [100], and [010]) corresponding to local min- 
ima of iii) the longitudinal mode with velocity 
Sli = \/ (cii -h Ci2 -|- 2c44)/2/9 polarized along [110] or 
along equivalent directions ([110], [110], and [110]) cor- 
responding to saddle points of "I>^^^ Accordingly, the 
oscillatory harmonics of the scattering rate are written 
as sums of three components: 
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The coefficients of order unity standing in the denomina- 
tors of these expressions are given by 
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where ca = 2c44-t-Ci2— Cn is the positive quantity charac- 
terizing the anisotropy of phonon spectrum. The contri- 
bution from transverse phonons is associated with piezo- 
electric potential, while the contribution from longitudi- 
nal phonons is due to deformation potential of electron- 
phonon interaction. 

The accuracy of the analytical approach leading to 
Eqs. (17)-(20) is not expected to be high for lowest- 
order PIRO peaks. Nevertheless, the applicability of 
these equations is steadily improved as one moves to 
lower magnetic fields and the argument of the oscillat- 
ing function becomes larger. For such low fields, the rate 
iy!^Q tends to overcome the rates and u^j^^. How- 
ever, in the samples with higher electron density, where 
the characteristic phonon wave number 2pF/fi is larger, 
the deformation-potential mechanism is more significant 
than the piezoelectric one. Therefore, both transverse- 
phonon and longitudinal-phonon contributions are im- 
portant and should be taken into account. Notice also 
that for low magnetic fields corresponding to the case of 
overlapping Landau levels one can take only the lowest 
harmonic A: = 1 in the sum in Eq. (13). This reduces the 
oscillating part of phonon-induced transport rate Vph to 
a simpler form 2d^j/^\ 

For completeness, one can also present the result ob- 
tained in the isotropic approximation, where there are a 
longitudinal mode and two degenerate transverse modes 
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with velocities si and St, respectively: 
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In this approximation, the phase of the oscillations, as 
well as the dependence on magnetic field and Fermi mo- 
mentum is different from the case of anisotropic phonon 
spectrum. 
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FIG. 1: (Color online) Phonon-induced resistivity calculated 
for the quantum well of experiment Ref. 17 at temperatures 
from 2 to 6 K. For T = 6 K, the partial oscillating contribu- 
tions from high-energy branch (LA) and low-energy branches 
(TA) are also shown (shifted for clarity). 

The numerical calculations described below use the 
following elastic constants for GaAs (in units 10^^ 
dyn/cm2): cn = 12.17 - 0.00144 r(K), ci2 = 5.46- 
0.00064 T(K), and C44 = 6.16-0.0007 r(K), taken from 
a semiconductor material database.^" The temperature 
dependence of these constants is not really essential in 
the range of T corresponding to experimental studies 
of oscillating magnetoresistance. The other GaAs pa- 
rameters used are T) = 7.17 eV, /114 — 1.2 V/nm, and 
p — 5.317 g/cm'^. An example of the results of numer- 
ical calculations according to Eq. (13) is presented in 
Fig. 1. The calculation employs the following param- 
eters of the high-mobility quantum well studied in Ref. 
17: ris — 3.75 X 10^^ cm~^, d^j = 30 nm, impurity-limited 
mobility 1.17 x lO'' cm^/V s, and quantum lifetime ow- 
ing to impurity scattering Tq — 15 ps. The theoretical 
plots are in agreement with experimental data of Ref. 17 
as concerns the background resistance, PIRO peak posi- 
tions, and the amplitude of the oscillations (for example, 
at T = 4 K the amplitudes of the peaks numbered 1 and 



2 are 0.16 and 0.075 Ohm per square, respectively). The 
calculation shows that the phonon-induced resistivity is 
determined by both longitudinal and transverse modes, 
whose contributions are nearly equal to each other at 
the fields S > 0.1 T and at 6 K. At lower fields the 
transverse- mode contribution becomes more important. 
With decreasing temperature both contributions are sup- 
pressed according to the factor TF{hu!\q/2T). The LA- 
mode contribution is suppressed stronger because the 
thermal activation energy for this mode is larger. Fig- 
ure 1 shows that the peak placed slightly below 0.6 T 
is the lowest-order one for the LA-mode contribution, so 
the origin of this peak (discussed in Ref. 17) seems to 
be clear. However, the relative amplitude of this peak 
is higher than that observed in the experiment, which 
is possibly related to overestimation of the deformation- 
potential interaction in comparison to the piezoelectric 
one.^'^ The other low-order peaks (n — 1, 2, 3) are formed 
by both transverse-mode and longitudinal-mode contri- 
butions, so their positions i?„ are not expected to follow 
exactly a l/S-periodic dependence. Nevertheless, the 
deviations appear to be small, and a 9-point linear fit of 
theoretical peak positions (see Fig. 2) gives the depen- 
dence n = 2pps/hijjc — 5n, where s ~ 3.44 km/s is very 
close both to the velocity sto — 3.40 km/s and to the ve- 
locity found experimentally,^'' while 5n ~ 0.21 is slightly 
larger than that obtained from the experiment. The 
parameter 5n determines the phase of the magnetoresis- 
tance oscillations, which appears to be close to the phase 
of the oscillating factor in the approximate analytical ex- 
pression (18), 5n — 1/4. Starting from T = 3 K, the 
theoretical dependence has a weak extra peak (marked 
by the arrow in Fig. 1) at B ~ 0.15 T, which comes from 
the longitudinal-mode contribution. Though the experi- 
mental magnetoresistance^^ does not show a peak in this 
region, a precursor of such a feature is apparently visible 
as a flattening of the minimum between the peaks n = 2 
and n = 3 at the temperatures above 4 K. 

It is important to compare the results of numerical cal- 
culation of the oscillating component of phonon-induced 
magnetoresistance with the results obtained from the an- 
alytical approach leading to Eqs. (17)-(20). Figure 3 
(a) shows the partial contribution of high-energy phonon 
mode (LA) calculated according to Eq. (13) together 
with the corresponding contribution based on the expres- 



sion I/. 



(fe) 

ph 



v"j^l , where h')^^ and v)^^ are given by 
Eqs. (19) and (20). In a similar way. Fig. 3 (b) shows 
both numerical and analytical, based on Eq. (18), re- 
sults for low-energy phonon modes (TA). The agreement 
between exact and approximate results becomes good in 
the low-field region corresponding to high-order magne- 
tophonon resonances. Both numerical and approximate 
results in Fig. 3 (a) demonstrate a beating pattern with 
a node at B ~ 0.08 T. The analytical consideration re- 
veals the origin of this beating, since it shows that the 
oscillations induced by LA phonons are formed as a sum 

(k) (k) 

of two contributions, i/j^q and Vj^i , with slightly different 
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FIG. 2: (Color online) Peak number n as a function of inverse 
peak position 1 / _B„ extracted from the calculated magnetore- 
sistance (Fig. 1) at 4 K. A linear fit for this dependence 
yields s ~ 3.44 km/s. The same kind of plot based on the 
experimental data is given in Ref. 17. The positions of low- 
order peaks slightly deviate from the linear dependence be- 
cause these peaks are formed by mixture of the oscillating 
contributions from different phonon modes. 



frequencies. 

Figure 4 presents the phonon-induced magnetoresis- 
tance calculated for the GaAs quantum wells with rig — 
10^^ cm~^ and = 13 nm studied in Ref. 14. Since the 
mobility of these samples (~ 6.6 x 10^ cm^/V s at T = 4.2 
K according to the resistance at i? = 0) is much lower 
than that in Ref. 17, both experimental^"* and theoreti- 
cal magnetoresistance show a few low-order oscillations. 
Because of large density and, consequently, large Fermi 
momentum in these samples, the main contribution 
to oscillating resistivity is caused by the deformation- 
potential interaction and comes from LA phonons. This 
is demonstrated by plotting partial contributions of dif- 
ferent modes in Fig. 4. At the temperatures 15-25 K all 
acoustic-phonon modes are fully activated, which means 
that the function F{huj\q/2T) is close to 1, so the re- 
sistivity should linearly increase with temperature. Ex- 
perimental plots of Ref. 14 indeed show a linear in- 
crease for the background resistivity, while the ampli- 
tudes of the oscillation peaks depend on T in a differ- 
ent way: the first peak is weakly modified by T and the 
second peak is suppressed with increasing T. Such a 
behavior can be explained by a decrease of the quan- 
tum lifetime with increasing temperature owing to in- 
elastic electron-electron scattering (see Ref. 17 and ref- 
erences therein). In the present formalism, this effect is 
described by adding the electron-electron scattering rate 
1/t|'^ — XT'^/hep, where A ^ 1, to the inverse quantum 
lifetime^^ I/t^ entering the density of states in Eq. (11). 
The result of the calculations improved in this way is 
also shown in Fig. 4. A reasonably good agreement with 
experimental data, as concerns the amplitudes of the os- 
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FIG. 3: (Color online) Partial oscillating contributions of the 
phonon-induced resistivity for the quantum well of experi- 
ment Ref. 17 at r = 4 K. (a) Contribution of high-energy 
branch (LA). Thick (black) lines correspond to numerical cal- 



culation of fp^' using Eq. (13), and thin (blue) lines are the 
results obtained from the approximate analytical expressions 
(19) and (20). The inset shows beating pattern in the low- 
field region, (b) Contribution of low-energy branches (TA). 
Thick (black) and thin (red) lines correspond to numerical 
calculation using Eq. (13) and to approximate analytical ex- 
pression (18), respectively. The low-field region is shown in 
the inset. 



dilations and their temperature dependence, is achieved 
at A = 2. 

The positions of the calculated magnetoresistance 
peaks deviate from the experimental data, though this 
deviation is not strong. The fitting procedure similar 
to that shown in Fig. 2 gives the characteristic velocity 
s — 5.1 km/s, which is smaller than s — 5.9 km/s de- 
termined in Ref. 14. It is worth noting, however, that 
the positions of the peaks in the experiment Ref. 14 are 
shifted to higher magnetic fields as the temperature in- 
creases. Though the reason of this shift is not clear, it 
signifies an increase of experimental velocity s with tem- 
perature. Therefore, one may expect a better agreement 
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FIG. 4: (Color online) Phonon-induced resistivity calculated 
for the quantum well of experiment Ref. 14 at different tem- 
peratures. For r = 15 K, the partial oscillating contributions 
from high-energy branch (LA) and low-energy branches (TA) 
are also shown (shifted for clarity). The additional (orange) 
plots take into account temperature-induced decrease of the 
quantum lifetime and are in a better agreement with experi- 
mental data of Ref. 14. 

of theory and experiment at lower temperatures. 

IV. SUMMARY AND CONCLUSIONS 

This paper presents a microscopic theory of magne- 
tophonon oscillations, also known as phonon-induced re- 
sistance oscillations (PIRO), in 2D electron systems. The 
phonon-induced resistivity pph is found by considering 
both deformation-potential and piezoelectric mechanisms 
of interaction of 2D electrons with bulk acoustic phonon 
modes. In the experimentally relevant situation of weak 
magnetic fields (large filling factors), the expression for 
the resistivity (or, eqivalently, for the phonon-assisted 
transport rate Vph) is determined by the electron-phonon 
coupling constants, phonon frequencies, and the density 
of electron states in magnetic field; see Eq. (13). The es- 
sential feature of the calculations is that the anisotropy 
of the phonon spectrum is explicitly taken into account. 
The theory is applied to quantum wells based on the ma- 
terials of cubic symmetry and grown in the [001] crystal- 
lographic direction. For this case, an analytical expres- 
sion for the oscillating part of the resistivity has been 



derived. This expression is valid in the region of fields 
corresponding to high-order magnetophonon resonances, 
when the ratio of characteristic phonon frequency to cy- 
clotron frequency is large. The calculations are carried 
out for GaAs quantum wells where PIRO have been ob- 
served experimentally. This includes the cases of 2D elec- 
tron gas with high density and moderate mobility^^ and 
with lower density and very high mobility. The calcu- 
lated magnetoresistance is in agreement with experimen- 
tal data. As concerns the positions of the magnetoresis- 
tance peaks, the agreement is very good in the case of 
high- mobility systems. 

Based on the results obtained, one may conclude 
that the bulk approximation for description of acoustic 
phonon modes and of their interaction with 2D electrons 
works reasonably good in application to magnetotrans- 
port in GaAs quantum wells. Another important con- 
clusion is that the velocity s in the empirical equation 
(1) does not, in general, correspond to a certain phonon 
mode. The reason for this is the complicated structure 
of the phonon spectrum. There are three anisotropic 
phonon modes (branches) with different velocities and 
two mechanisms, deformation-potential and piezoelec- 
tric, of electron-phonon interaction. Relative contribu- 
tions of these modes and of the interaction mechanisms 
depend on 2D electron density, temperature, and other 
parameters such as the width of the quantum well and 
growth direction. This explains why the velocity s deter- 
mined by fitting the experimental PIRO peak positions 
is expected to vary in different experiments, especially 
when such a fit is based upon a few low-order peaks 
which have the highest amplitudes and, therefore, are 
best visible experimentally. Indeed, the low-order peaks 
are typically formed as a superposition of contributions 
from different modes, and even for each single mode there 
is a mixture of phonons with different velocities owing to 
the anisotropy. The anisotropy influences not only the 
frequency but also the phase and the amplitude of the 
oscillations. 

On the other hand, in the region of high-order mag- 
netophonon resonances the resistance oscillations are de- 
scribed by a fixed set of phonon velocities determined 
here for the case of [001]-grown quantum wells; see Eqs. 
(17)-(20). The phases of these oscillations also become 
definite. Moreover, if the density of 2D electrons is not 
large, so the piezoelectric-potential interaction is signif- 
icant, the main contribution to magnetoresistance oscil- 
lations in [001]-grown wells is characterized by a single 
velocity s = sto = \/c^aIp corresponding to a TA mode 
propagating in the quantum well plane and polarized per- 
pendicular to this plane. This is the case of experiment 
Ref. 17, where many magnetoresistance oscillations have 
been observed owing to a very high electron mobility. 
For the wells whose growth axes are different from [001] 
the characteristic phonon velocities have to be different. 
These cases are not studied here in detail because PIRO 
have not yet been observed in such wells. Nevertheless, 
based on the consideration presented in Sec. Ill, the fol- 
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lowing general method for determination of characteristic 
phonon velocities in the wells of arbitrary growth axis is 
proposed. Take the surfaces of equal frequency w defined 
for each mode in the q-spacc by the equation sl^^q = oj. 
Find the points (numbered by the index i) where these 
surfaces touch the surface of the cylinder of radius 2pF/h 
and axis along the growth direction (this cylindrical sur- 
face is given by the equation qsinx = 'i.pp/K). The 
set of frequencies u = u)\i (or, equivalently, velocities 
s\i = cjxi/q) corresponding to these points will give a 
set of oscillating harmonics with arguments 27Tkoj\i/oJc 
in the magnetoresistance. An equivalent procedure is the 
search for extrema of the function / sin x, as described 
in Sec. III. To find which frequencies out of oj\i give the 
main contribution and to determine the phases and am- 
plitudes of the oscillations, a detailed analysis is neces- 
sary. In the case of [001]-grown wells a single TA mode 
gives the main contribution [Eq. (18)] because the sur- 
face of equal frequency for this particular mode touches 
the cylinder not just in a finite number of points, but 
along the whole circumference of the cylinder at x = 7r/2. 
To study the influence of different phonon modes on 



the magnetoresistance oscillations, it is desirable to carry 
out experiments in high-mobility 2D systems with den- 
sities varied in the range 4 — 10 x 10^^ cm~^. If it is 
possible to reach high mobilities for the structures with 
growth axes different from [001], measurements of the os- 
cillating magnetoresistance in such systems would be also 
important for investigation of the phonon anisotropy ef- 
fects discussed above. The author hopes that the present 
theoretical work may stimulate such experiments. 

The consideration given in this paper assiimes the lin- 
ear response regime, when both electron and phonon sys- 
tems arc; close to eqiiilibrium. One may expcict a more 
interesting and rich behavior of the phonon-induced resis- 
tance oscillations under nonlinear transport regime when 
the electric current through the sample increases. Studies 
in this direction are already undertaken. ^^'^^ In addition 
to the work already done, more efforts, both experimental 
and theoretical, are necessary for better understanding of 
the mechanisms of phonon-induced oscillations and their 
interplay with the other microscopic mechanisms respon- 
sible for the magnetoresistance of 2D electrons at large 
flUing factors. 
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